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THE DISTRIBUTION OF ODD NITROGEN IN THE LOWER 
STRATOSPHERE AND POSSIBLE PERTURBATIONS CAUSED BY 
STRATOSPHERIC AIR TRANSPORT 


Ivar S. A. Isaksen and Eigil Hesstvedt 
Institute of Geophysics 
University of Oslo 
Oslo, Norway 


Abstract 


In the lower stratosphere a significant production of odd nitrogen results 
from the reaction N 2 0 + O ( *D) -> 2 NO. Since the transport is relatively slow, 
odd nitrogen builds up with a maximum mixing ratio of 2 x 10 * 1 * * * * * * 8 at 30km. Pro- 
files of odd nitrogen, for different latitudes, winter and summer, are computed 
from one-dimensional transport models. Variations with latitude are small. 
Horizontal transport is therefore not believed to alter our results significantly. 

In order to evaluate the effect of odd nitrogen upon the ozone layer, NO x 
profiles are calculated. OH is here a key component, since it converts NO z to 
HNO 3 . In the region where ozone is determined by chemistry rather than by 
transport (above 25km), N0 2 is found to be relatively abundant. 

The effect of stratospheric transport on the NO x distribution is shown to 
depend critically upon the height of emission. The effect increases by a factor 
of 5 or more for a change of flight level from 18 km to 23 km. This strong de- 
pendence should be duely considered when future stratospheric transport is 
discussed. 


1. Introduction 

It is well established (C rut zen, 1971, Johnston, 1971, Hesstvedt, 1973), 

that nitrogen oxides act upon the ozone layer through the reaction 

N0 2 + O ( 3 P) -► NO + 0 2 kj = 9.2x 10' 12 (1) 

whereby odd oxygen is destroyed. Accordingly, reliable information about the 

distribution of nitrogen oxides in the lower stratosphere is essential in quantita- 

tive assessments of this effect. Production of odd nitrogen takes place in two 

regions (Crutzen, 1971, Isaksen, 1973). In the upper region, around 100km, 
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ionic processes result in an important. production of NO and N. In the lower re- 
gion, between 20 and 30km, NO is produced when N 2 O is destroyed by O ( 1 D). 

In this region the chemistry appears to be rather complicated. Modeling of the 
atmosphere with respect to odd nitrogen species, must also comprise a series 
of oxygen, hydrogen, and carbon species. The model has earlier been presented 
by Isaksen (1973) and will therefore not be disucssed in detail here. We will, 
however, discuss production and loss reactions with O ('D) and OH, which are 
day components for the nitrogen oxides in the lower stratosphere. 


2. The model 

Important parameters in the NO x calculations, such as eddy diffusion coef- 
ficient, temperature and the air density, are known to vary markedly with sea- 
son and latitude. Calculations are therefore made for 60° and 30°, winter, for 
the Equator, and for 30°, 60°, and 90°, summer. Temperatures and densities 
are taken from CIRA (1965). The models are one-dimensional, and vertical 
eddy transport is included. The eddy diffusion coefficients used are given in 
Figure 1. Up to 27km the coefficients are taken from Gudiksen et al. , (1968). 
Above 27km an exponential increase is assumed, as suggested by Lindzen(l971). 
His indication of a seasonal variation is also adapted. 

Hesstvedt (1973) has shown that a two-dimensional model should be con- 
sidered for computations of the ozone distribution in the lower stratosphere. 
Since our model is one-dimensional we therefore decided to use observed values 
of ozone (Hering and Borden, 1964). 

Our model is a steady state model, where the photochemistry is averaged 
over the day. Loug lived species which are influenced by eddy diffusion are de- 
termined from equations of the type 


d [X] 
dt 


d 

dz 


r n d5 x , 

K z [M] — ) + P, 
dz 


Qx [x] = 0 


( 2 ) 


where K z is the vertical eddy diffusion coefficient, M the total number density, 

8 X - [X] / [M] the mixing ratio of component X, and P x and Q x [X] are averaged 
photochemical production and loss terms . This equation is solved by the direct 
implicit method given by Isaksen (1973), using finite differences. 

Components with short chemical lifetimes, hours or less, are not influenced 
by transport processes. The transport term in equation (2) may then be disre- 
garded, and average daytime values are determined from expressions of the 
type 

Px ~ Qx [X] = 0 (3) 

where P x and Q x [x] are average daytime production and loss terms . 
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3. Atomic oxygen in the excited state and hydroxyl 


Two components play an important role in the formation of odd nitrogen 
species; these components are O ( x D) and OH. Odd nitrogen is produced 
through the reaction ( 1 ) which appears to be the only effective source in the 
stratosphere. Nitrogen oxides will therefore depend strongly on the O ( x D) pro- 
files , which are determined from 

0 3 + hi> -► 0 2 + O ( 1 D) X < 311 nm (4) 

O ( J D) + M ->0 ( 3 P) + M k 5 = 6 x 10 - 11 (5) 

which gives 

[0( 3 D)] = J 3 • [0 3 ] A 5 • [M] ( 6 ) 

In our model observed values are used for O 3 and M. J 3 is the daytime average 
dissociation rate of 0 3 for X < 311 nm. Its strong variation with height in the 
lower stratosphere (Isaksen, 1973) is a result of ozone absorption, and can 
therefore be calculated from the observed ozone profile. 

O (*D) profiles, calculated from equation ( 6 ) are given in Figure 2 for dif- 
ferent latitudes, summer and winter. The rapid decrease.below 50km is a re- 
sult of the decrease in J 3 . Variations with latitude are small except at high 
latitudes, winter, where the production is very slow because of the sun. 

Hydroxyl is a second component in the lower stratosphere. It is produced 
by the reaction 

. Moo 

NO + H0 2 -* N0 2 + OH k 7 = 5 X 10 ' 11 • e' T (7) 

In addition to this reaction there is some production through the reaction 

O ( 3 P) + H0 2 -* OH + 0 2 k 8 = 7 x 10 - 11 ( 8 ) 

Below 30km, however, this reaction becomes negligible, since O ( 3 P) de- 
creases rapidly below this height. OH is lost through several reactions. Im- 
portant reactions in the lower stratosphere are found to be 

_ 1886 

OH + CH 4 -*H 2 0 + CH 3 k 9 = 5. 5 x 10 -13 • e T (9) 
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_ 115 

OH + CO -*H + C0 2 k 10 = 2.1 x 10“ 13 • e T (10) 

. 1650 

OH + HN0 3 H 2 0 + N0 3 kn= 1.7 xlO" 11 • e T (11) 

. 170 

OH + N0 2 + M ->-HN0 3 + M k 12 = 1.05 x 10' 11 • e T (1 + f)' 1 

(12) 

f = 1. 6 x 10’ 2 • T 4 / (4 x lO" 11 [M] ) 

Reactions (9) and (10) are real loss reactions, since CH 3 and H are ultimately 
converted to H0 2 . Reaction (12), on the other hand, is not a real loss since 
most of OH lost by this reaction is reformed through the reaction 

HN0 3 + hi/-*OH + N0 2 X < 546nm, J aHN03 (13) 

The real OH loss is only the fractional part going through reaction (11) and 
reaction 


HN0 3 + h i> -»■ HO 2 + NO X< 370nm, JbHN0 3 (14) 


on the short wavelength side of HN0 3 dissociation. It appears that reaction 
(14) is much less effective than reaction (13) . OH may then be determined 
from 


[OH] = 


(k 7 [no] + k 8 [o ( 3 P) ] ) [ho 2 ] 
k 9 [CH 4 ] + K 10 [CO] + 2k n [HN0 3 ] + J bHN Q 3 


(15) 


OH can therefore only be determined in a model where all the components oc- 
curring in equation (15) are computed simultaneously. This is done in our 
model. Calculations are repeated until convergence is obtained. The resulting 
OH profiles are shown in Figure 3. OH has, for all latitudes, a maximum near 
45km with a peak value of 2 - 3 x 10 17 cm‘ 3 . Down to the peak level H0 2 is ef- 
fectively converted to OH through reaction (8) , while the destruction through 
reactions (9), (10), and (12) is slow. Below 45km this is no longer true and we 
get a rapid decrease in OH down to the tropopause, where we find a minimum 
value near 10 s cm -3 . The strong variations with latitude below about 17km 
come from the variations in the tropopause level. At the Equator heights below 
17km are in the troposphere, at high latitudes all heights above 8km are in the 
stratosphere . 
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4. Distribution of nitrous oxide 


Destruction of N 2 0 through the reaction 

N 2 0 + O (*D) -> NO + NO k 16 = 9 x 10" 11 (16) 

is the only important source of NO x in the lower stratosphere. Realistic infor- 
mation on the spatial distribution of this component is therefore essential. N 2 0 
has been found to be fairly uniformly distributed in the troposphere with negli- 
gible seasonal variations. A mixing ratio of 2. 5 10 ~ 7 is regarded as represen- 
tative for the troposphere (Schiitz et al. , 1970). Production of N 2 O is expected 
to be through soil bacteria. So far no photochemical sources are known. It is 
therefore difficult to estimate production rates. Since the loss rate is small, 
vertical transport maintain a constant mixing ratio of N 2 O in the troposphere, 
and the above value is used as boundary value at the tropopause. In the lower 
stratosphere the conditions are quite different. Diffusion is much lower and 
loss becomes more effective through the reaction 

N 2 0 + hi>->N 2 + O X < 337. 5nm (17) 

and through the slower reaction 

N 2 0 + O ( : D) ^N 2 + 0 2 k is = 9 x 10' 11 (18) 

in addition to reaction (16). N 2 0 profiles are therefore obtained from 

■^( Kz [M] Tt. (®h 2 o)) - ( j n 2 o +(ki6 + k,8) [o(‘D)]j [n 2 0] =0 (19) 

The results are shown in Figure 4. The slope of the profiles depends on the 
dissociation rates and on the eddy diffusion coefficient. Weak diffusion leads to 
a rapid decrease in the N 2 O mixing ratio (5n 2 o) . This occurs at 25-30km, 
where K z has a minimum. Variations with latitude and season are relatively 
small, except for high latitude, winter. Apart from this case, production of 
NO x through destruction of N 2 O is confined to levels below 3 5-40 km. At 60° , 
winter, the low sun makes photodissociation less important, and when in addi- 
tion the eddy mixing is stronger, N 2 0 is transported to greater heights before 
it is dissociated. 
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5. Production and distribution of odd nitrogen 


From the profiles of N 2 0 and O ( ! D), shown in Figures 2 and 4, the produc- 
tion rates of odd nitrogen are easily calculated from 

P NOy = 2k 16 • [0(»D)] [N 2 0] (20) 

The efficiency of this production depends on how efficient reaction (16) is in 
breaking down N 2 0, compared to reactions (17) and (18). Figure 5 shows the 
loss rate of N 2 0 from reactions (16) and (17). Loss by reaction (18) is equal to 
the loss by reaction (16). The curves are for 30° summer. It is clearly demon- 
strated that reaction (17) is the main loss reaction of N 2 0. Only a few percent 
of the loss will produce nitrogen oxides through reaction (16). 

Nitrogen oxide production from eq. (20) is given in Figure 6. In the summer 
production varies relatively little with latitude. Differences below 20km is a 
result of the variation in tropopause level, which influences the O (*D) densities. 
At all latitudes , odd nitrogen production increases rapidly up to a marked max- 
imum between 22-26 km. In the winter the production is important over a much 
wider region. This is clearly demonstrated for 60°, winter, where a substan- 
tial production takes place above 40km, with a maximum around 35km. The 
winter profiles are a result of the combined effects of the absorption of radia- 
tion producing O ( J D) and the strong eddy diffusion, which effectively transports 
N 2 0 to greater heights. At 30°, winter, these effects are less pronounced, with 
the result that we get a profile more like the summer profiles . Maximum pro- 
duction rates vary from 1. 5 x 10 2 cm -3 s -1 in the summer hemisphere, to 1. 5 x 
10 1 cm' 3 s‘' at 60° , winter. If we add up the total production in a column, the 
average production for the whole globe becomes 10 8 cm _2 s _1 . 

The distribution of total odd nitrogen in the lower stratosphere is easily 
calculated when its production rate is known. The only photochemical loss of 
nitrogen oxides in the atmosphere is through the reaction 

NO + N^N 2 + 0 k = 2.2 x 10 -11 (21) 

However, atomic nitrogen decreases rapidly below 60km (Isaksen, 1973). This 
reaction is. thereto re of no interest for our problem. Odd nitrogen is therefore 
removed from the lower stratosphere by vertical eddy diffusion alone. Since 
K z increases with height, the particle flux of odd nitrogen, NO y 

F = -K z [M] £ (s NOy ) (22) 

is directed downwards. 
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The flux is proportional to the eddy diffusion coefficient and to [m] . It is, 
therefore clear that the increase in 5 NOy , sufficient to transport away the odd 
nitrogen produced at a certain height, depends strongly on the heights where the 
production takes place. Assuming steady state conditions, 



gives the height distribution of nitrogen in all odd forms, i.e. , in the forms of 
NO, NO2 , and HNO3; all other odd components are present in negligible quan- 
tities. Results of this calculation are given in Figure 7. As might be expected 
from the considerations above, there are latitudinal variations in odd nitrogen, 
due to variations in K z and in production rates. Above 20km, where the effect 
on ozone is of importance, variations with latitude are very small in the sum- 
mer. The mixing ratio of odd nitrogen increases strongly up to 25km, the re- 
gion of the strongest production (Fig. 6). But even above this height there is 
substantial increase to above 30km, partly due to the low K z values. Above 
35km Sjvjoy drops as a result of the low mesosphere values used for NO(Isaksen, 
1973 a). Maximum mixing ratios of 1.2 x 10' 8 are found around 32km. 

In the winter the low solar elevation leads to production at higher altitudes . 
Therefore the increase in mixing ratio takes place at higher altitudes. This is 
especially true for 60° , winter, where the strong atmospheric absorption gives 
high mixing ratios in the upper stratosphere. It is, therefore, an increase in 
odd nitrogen mixing ratios in the whole stratosphere. 


6. The distribution of odd nitrogen on NO, N0 2 , and HN0 3 


In order to estimate the effect on ozone, we shall next see how the total odd 
oxygen, as computed above, is distributed on the different odd nitrogen species. 
The initial product is NO through reaction (16). This component will rapidly be 
converted to NO 2 through reaction 

( I no r \ 

— J (24) 

The balance between NO and N0 2 is maintained through solar dissociation of 

no 2 


N0 2 + hi> -► NO + O ( 3 P) 397. 5nm (25) 

with some contribution from reaction (1). 
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The conversion to HN0 3 is much slower (reaction (12)). To obtain HNO 3 
densities, the loss reactions must be compared with reaction (12). The main 
loss of HNO 3 is through reactions (13) and (14) with some contribution from re- 
action ( 11 ). N0 3 is rapidly converted to NO and N0 2 through the reactions 


NO 3 

+ hy -*■ NO + 0 2 

X < 11/i 

NO 3 

+ hi> -* N0 2 + 0 

X < 700nm 


(26) 


NO 3 is unlikely to be of importance and will be disregarded here. 


From the above reactions the distribution of odd nitrogen on NO, NO 2 , and 
HNO 3 can be calculated. Equilibrium between NO and NO 2 is given by reactions 
(24), (25), and ( 1 ). Loss rates of NO and NO 2 are shown in Figure 8 for 30° 
summer. It is seen that equilibrium is established after a few minutes. Sim- 
ilarly, equilibrium between NO - N0 2 and HNO 3 is given by equations ( 11 ) - (14). 
Equilibrium in this cycle is much slower, as shown in the same figure. Life- 
times decrease from 1 day at 30 km to 10 days at 15km. The short lifetimes in 
the NO to NO 2 cycle means that NO is converted to NO 2 a few minutes after sun- 
set, when there is no photodissociation. Likewise, NO is also reproduced 
rapidly after sunrise from photodissociation of N0 2 . 

The long lifetimes in the NO - NO 2 - HNO 3 cycle result in negligible HNO 3 
variations over the day below 30km. By confining the equations for equilibrium 
in the NO - N0 2 and in the NO - N0 2 -HN0 3 cycle, average daytime number 
densities for NO, N0 2 , and HN0 3 are given by 


[NO] = [ NOy ] 


Jhno 3 + kn [OH] 

Jhno 3 + kn [OH] +k ' 12 [OH] [M] 

k i [0( 3 P)] +J n02 

k l [°] + J N0 2 + k 24 [0 3 ] 


(27) 


^hno 3 + k ll t OH ] 

[NO 2 ] = [ NO ] 

j hno 3 + k n [OH] + k ' 12 [OH] [M] 

k 24 [O3] 


k i [0( 3 P)] + J no 2 + k 24 [o 3 ] 
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( 29 ) 


k 12 [OH] [M] 

[HN0 3 ] = [.'NO ] 

Jhno 3 + k u [OH] + k ' 12 [OH] [M] 

k [ 2 gives the efficiency of the coefficient k 12 in converting NO + N0 2 to HN0 3 : 

, k 24 [° 3] 

k 12 =k 12 (30) 

k 24 [o 3 ] +kj. [0(3p)] + j N02 

The relative abundances of the three nitrogen species, given by equations 
(27), (28), and (29), are presented in Figure 9 for heights between 10 km and 
40km for 30° summer. HN0 3 is seen to be the dominating component below 
about 27km. At 25km 70% of odd nitrogen present is in the form of HN0 3 , and 
it exceeds 40% at all heights below. The maximum HNO 3 concentrations around 
25km is a result of the high ozone concentration. Below about 15km HN0 3 in- 
creases due to the increase in OH. Above 35km HN0 3 is present in small quan- 
tities. N0 2 , which is of special interest in the ozone chemistry, is most abun- 
dant around 30km, where 50% of the total odd nitrogen is in the form of N0 2 . 
Relative densities vary little below 25km, between 20% and 30%. NO is the 
most abundant component above 35km, compared with the other components. 
Below 30km it is a minor component. The NO - N0 2 cycle is of interest for 
ozone destruction through reaction ( 1 ). As we have already seen, lifetimes are 
very short. During the day, when ozone is lost through reaction (1), NO and 
N0 2 are in photochemical equilibrium. When the sum of NO and N0 2 is known, 
equilibrium values are easily calculated from equations (1), (24), and (25). 

NO x mixing ratios are given in Figure 10 for summer and winter hemi- 
sphere. Above 30 km the mixing ratios differ little from odd nitrogen mixing 
ratios, due to low HN0 3 densities (Fig. 8 ). Below 25km, there is a sharp drop 
down to 10 km, where mixing ratios of about 10" 7 * 9 are obtained for all latitudes. 
Variations with latitude are very small, except in the winter hemisphere where 
S NOx falls off at higher latitude. 

7. Air transport in the stratosphere and its effect on the NO x distribution 

It is evident that nitrogen oxides produced from aircraft flying in the strato- 

sphere will alter the stratospheric NO x distribution. Since vertical eddy diffusion 
varies strongly with height, it is of interest to estimate the magnitude of this 
additional source when the aircraft are assumed to fly at given heights. We shall 
here consider the effect of a "reference fleet" of 200 four-engine aircraft, each 
performing 1000 flights of 2 - 1/2 hours per year. The fuel flow is assumed 
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to be 5000kg per hour per engine. The emission of NO (converted to N0 2 ) is 
assumed to be 21 g per kg fuel. These 200 aircraft are assumed to fly at the 
same flight level uniformly routed around the globe. This gives an average pro- 
duction rate of 1. 88 x 10 7 NO x molecules cm _2 s _1 . Since we use finite differ- 
ences with steps of 1km, we have, in our model, a production rate of 1.88 x 
10 2 NO x molecules cm- 2 s _1 . This number will be used as an additional source 
in equation (23) at given flight levels. The variation in <5 NOx , for flight levels 
between 15km and 32 km, is given in Figure 11. The steady state profiles in 
Figure 10 for 60°, summer, and 30°, winter, are used as a basis, and the in- 
crease is given in percent. Since the ozone distribution below 25km is largely 
dominated by transport processes (Hesstvedt, 1973 a), the increase in NO x be- 
tween 25km and 30km is of primary interest for the ozone problem. 

As might be expected, the increase in NO x depends strongly on the flight 
level. NO x injection at high altitudes leads to a considerable increase in NO x . 
This can be explained by the very fast drop in eddy diffusion coefficient above 
the tropopause (Fig. 1). The added NO x has to be transported downwards into 
the troposphere. Production at heights where K z is small, therefore means that 
the gradient in Sno x must be much greater than the gradient is at lower heights 
with higher K z values, in order to maintain the flux of NO x given by equation 
( 22 ). 

The increase in 5 no x is somewhat different in the two diagrams of Figure 
11, a result of difference in K z . In both cases, however, the perturbation in 
NO x increases rapidly as the source level, given along the abcissa, is moved 
upwards. With the source at 18km, the increase in NO x is approximately 2% 
for 60°, summer, and 0. 5% for 30° , winter. With the source at 23km, NO x in- 
creases with 5-10% in the 25-30km region in both models, and with source at 
28km, NO x increases by around 20% in both models. The calculations give sim- 
ilar results for the other latitudes. These results have been used in a two- 
dimensional model to show the effect on the ozone layer (Hesstvedt, 1973 a). 

These numbers clearly demonstrate that the effect of stratospheric trans- 
port on the NO x densities depends strongly on the flight levels. Flight levels 
above 20km lead to a considerable increase in NO x . Having in mind that the 
perturbation (in %) of the ozone layer is roughly one third of the perturbation 
(in %) of NO x , the importance of the flight level should be duely considered when 
the effect of stratospheric transport upon the ozone layer is discussed. 
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Vertical eddy diffusion coefficient. 

Vertical profiles of O (*D) for 90° summer, 60° summer, 30° 
summer, 0° , 30° winter, and 60° winter. 

Vertical profiles of OH, for the same latitudes as in Figure 2. 

Vertical profiles of N 2 O for 90° summer, 0°, 30° winter, and 
60° winter. 

Vertical profiles of N 2 0 loss rates from equations (16) and (18) 
at 30° summer. 

Vertical profiles of production rates of odd nitrogen from 
equation (16) for the same latitudes as in Figure 4. 

Profiles of odd nitrogen mixing ratios for the same latitudes 
as in Figure 2. 

Vertical profiles of odd nitrogen loss rates in the NO - N0 2 
cycle (curves q N02 “ no> 311(1 3 no) 311(1 in the NO - N0 2 - HNO 3 
cycle (curves qN0 2 - hno 3 311(1 <1hno 3 ) for 30° summer. 

Distribution of odd nitrogen on NO, N0 2 , and HN0 3 . 

Vertical profiles of NO x mixing ratio. 

Increase (in %) of NO x mixing ratio caused by a fleet of 200 
SST aircraft as a function of flight level (horizontal axis). 
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Vertical eddy diffusion coefficient, cm 

Figure 1. Vertical eddy diffusion coefficient. 
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Number density (particles cm' 3 ) 

Figure 2. Vertical profiles of 0 1 1 D) for 90° summer, 60° summer, 30° summer, 0°, 30° winter, and 60° winter. 
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Number density ( particles cm" 3 ) 

Figure 3. Vertical profiles of OH , for the same latitudes as in Figure 2. 
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Figure 5. Vertical profiles of N 2 O loss rates from equations (16) and (18) at 30° summer. 
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Mixing ratio 

Figure 7. Profiles of odd nitrogen mixing ratios for the same latitudes as in Figure 2. 
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Figure 9. Distribution of odd nitrogen on NO, NO 2 , and HN0 3 . 
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N0 X mixing ratio 

Figure 10. Vertical profiles of NO x mixing ratio. 
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Flight level, km 

Figure 1 1. Increase (in %) of N0 X mixing ratio caused by a fleet of 200 SST aircraft as a function of flight level (horizontal axis). 
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